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Abstract 


The Advanced Automated Directional Solidification Furnace (AADSF) is a Bridgman-Stockbarger micrograv- 
ity processing facility, designed and manifested to first fly aboard the second United States Microgravity Pay- 
load (USMP-2) Space Shuttle mission. The AADSF was principally designed to produce high axial thermal 
gradients, and is particularly suitable for metals solidification experiments, including non-dilute alloys. To ac- 
commodate a wider range of experimental conditions, the AADSF is equipped with a reconfigurable gradient 
zone. The overall design of the AADSF and the relationship between gradient zone design and furnace per- 
formance are described. Parametric thermal analysis was performed and used to select gradient zone design 
features that fulfill the high thermal gradient requirements of the USMP-2 experiment. The thermal model and 
analytical procedure, and parametric results leading to the first flight gradient zone configuration, are pres- 
ented. Performance for the USMP-2 flight experiment is also predicted, and analysis results are compared to 

test data. 


Introduction 

The NASA Marshall Space Flight Center con- 
ceived and funded development for the Advanced Auto- 
mated Directional Solidification Furnace (AADSF), 
which will fly on the second United States Microgravity 
Payload (USMP-2) mission in early 1994. The AADSF 
is a Bridgman-Stockbarger crystal growth furnace de- 
signed to operate in a microgravity environment on or- 
bit. Manifested on several future USMP science 
missions, the AADSF flies on a multi-purpose experi- 
ment support structure (MPESS) in the space shuttle 
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Figure 1: AADSF System Flight Configuration 

•Branch Manager, Space Programs Division 
Member A1AA 


payload bay. For the first flight experiment, the AADSF 
is configured to maximize the axial temperature gradi- 
ent at the crystal growth interface. This paper describes 
the engineering process used to select the gradient zone 
design features that enable the AADSF to produce the 
high thermal gradients needed by the USMP-2 
experiment. 

The AADSF system consists of three major 
components, the furnace assembly, the Data Acquisition 
System computer (DAS), and the Signal Conditioning 
and Control System (SCCS). Figure 1 depicts the 
AADSF flight elements mounted on the MPESS. A 
sealed experiment apparatus container (EAC) provides 
an inert gas atmosphere for the three zone Bridgman- 
Stockbarger furnace and the sample translation and ex- 
periment telemetry subsystems. The DAS collects ana- 
log data for downlink during processing. The SCCS 
conditions sample and control thermocouple signals, 
stores experiment timelines, and accomplishes closed- 
loop control of the individual furnace heaters and the 
sample translation system. 

Furnace Design 

The AADSF can ultimately attain hot zone 
temperatures of up to 1600°C. However, for the first 
two flights, heat extraction plate and control thermo- 
couple material choices limit the maximum operational 
temperature to 1 150°C in the hot zone and 850°C in the 
cold zone. The furnace bore diameter permits process- 
ing with sample diameters up to 20 mm. The furnace 
includes a 250 mm long isothermal hot zone, and a 100 
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mm cold zone. The system accommodates sample 
translations up to 240 mm. The ability to change 
insulation and heat extraction plate thickness and mate- 
rial makes the gradient zone reconfigurable {1]. Chang- 
ing the gradient zone adapts the furnace to different 
experimenter requirements. Figure 2 shows a section 
view of the furnace. 



Figure 2: AADSF Cross Sectional View 


To produce high temperature gradients within 
a sample, heat transfer across the gradient zone into the 
cold zone is minimized. The heat extraction plate ac- 
complishes this by conducting heat from the hot zone to 
the furnace cooling system before it can reach the cold 
zone. Low thermal conductivity gradient zone insula- 
tion further isolates the cold zone from the hot zone. 
For the USMP-2 mission, the engineering performance 
specification required that the furnace produce tempera- 
ture gradients of 100°C/cm in an instrumented quartz 
rod with the hot zone operating at 1 100°C. This design 
criterion ensured sufficient temperature and gradient 
performance for the planned first flight experiment. 


with anticipated hot and cold zone temperatures near 
850°C and 350°C, respectively. 

Unlike a conventional two-zone Bridgman fur- 
nace, the AADSF has a booster heater at the interface 
between the hot zone and the gradient zone. This inde- 
pendently controlled heater maintains hot zone iso- 
thermality near the gradient zone. Without the booster 
heater, heat transfer from the hot zone to the heat ex- 
traction plate would cause furnace bore temperatures to 
decrease near the gradient zone, thereby degrading iso- 
thermality. The booster heater reinforces the hot main 
heater in this local region, improving hot zone iso- 
thermality and increasing the axial thermal gradients 
produced within the sample. 

Crystal growth in the AADSF occurs with the 
translation of the sample from the hot end of the fur- 
nace through the gradient zone, where directional so- 
lidification occurs. The translation system pulls the 
sample through the gradient zone at controlled rates be- 
tween 0.5 and 50 mm/hour. The translation rate is se- 
lected by the experimenter, based on the expected 
crystal growth parameters, and may be changed during 
processing without significantly disturbing the thermal 
performance of the furnace. 

Experiment Design 

The first flight experiment using the AADSF 
will be "Growth of Solid Solution Crystals - HgCdTe" 
aboard the USMP-2 Space Shuttle flight. The principal 
investigator for this experiment is Dr. S. L. Lehoczky, 
with co-investigators Dr. F. R. Szofran and Dr. D. C. 
Gillies, all of the NASA Marshall Space Flight Center. 

For the USMP-2 experiment the AADSF sys- 
tem will process directional solidification of Hg^Cd/Te 
crystals under microgravity and high thermal gradient 
conditions. This material behaves as a nondilute alloy 
of HgTe and CdTe. Under normal earth gravity condi- 
tions, this material resists compositionally uniform 
growth by the Bridgman-Stockbarger method. As the 
material solidifies, CdTe is preferentially incorporated 
into the crystal, and the liquid sample becomes enriched 
with denser HgTe near the solidification interface. Un- 
der normal earth gravity conditions, the vertical Bridg- 
man method results in increased HgTe enrichment near 
the crystal/melt interface, because nondilute alloy com- 
ponents tend to separate due to density differences. Un- 
der microgravity conditions, diffusion becomes the only 
significant mechanism for HgTe enrichment in the 
melt, and density stratification should be nearly nonex- 
istent, hence improving the compositional uniformity of 
the crystal. 

The nondilute nature of Hg^CdJe also drives 
the high gradient requirement of the USMP-2 
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experiment. High temperature gradients within the 
sample result in better control of the ciystal-melt inter- 
face. High axial gradients tend to flatten the shape of 
the solidification interface. Also, constitutional super- 
cooling, a problem that can result in polycrystalline 
directional solidification of nondilute alloys, is avoided 
by sufficiently high thermal gradients coupled with suf- 
ficiently low growth rates. Compositional variations in 
the liquid sample also cause variations in the solidifica- 
tion temperature during the crystal growth process. 
Such solidification temperature variations cause axial 
displacement of the interface, and high temperature 
gradients minimize this displacement [2]. 

Gradient Zone Design Options 

Design options for the AADSF gradient zone 
included: (1) gradient zone material and thickness, and 
(2) heat extraction plate material and thickness. Table 1 
lists properties for the candidate gradient zone materi- 
als. The USMP-2 experimenter preferred the shortest 
possible gradient zone, commensurate with achieving 
hot and cold zone temperature control at 850°C and 
350°C, respectively. Based on this requirement and ear- 
ly analysis, zone insulation thicknesses of 7.6 mm and 
4.4 mm were selected to separate the heat extraction 
plate from the hot zone and cold zone heater cores, re- 
spectively. In addition to the insulation components, the 
gradient zone also includes a 0.5 mm thick washer 
separating the cold zone heater core and gradient zone 


Material 

Density 

Thermal 
Conduc- 
tivity (at 
1000 °C) 

Specific 
Heat (at 
1000°Q 

Maximum 

Use 

Temperature 


(g/cm 1 ) 

(W/m-°C) 

(J/kg-«C) 

(°Q 

Gradient Zone Insulation Material Candidates 

FBD™ 

fibrous 

1.4 

033 (ave) 

657 

2^200 

zirconia 





AL-30™ 

fibrous 

alumina 

0.5 

0.2 

1,050 

1,540 

ZAL-45™ 

fibrous 

alumina 

0.7 

03 

1,050 

1,650 

Alumina 

(99.5% 

dense) 

3.9 

6.0 

1,230 

2,050 

Heat Extraction Plate Material Candidates 

Platinum 

21.5 

84 

159 

1,600 

Inoonel718 

8.2 

27 

712 

1,100 


(FBD. AL-30. and ZAL-45 are trademarks of Zircar Products Inc.) 

Table 1 - Properties of Gradient Zone Insulation and Heat 
Extraction Plate Material Candidates 


insulation, and the heat extraction plate. Figure 3 de- 
picts the details of the gradient zone design. 



Figure 3: AADSF Gradient Zone Detail 

With the thicknesses defined, gradient zone in- 
sulation material remained to be selected. Candidate 
gradient zone insulation materials included several low 
density ceramic fiber insulations and dense alumina. 
Pure alumina and alumina/zirconia blends were the 
principal candidate fibrous insulation materials. Poten- 
tial silica contamination of platinum alloy control ther- 
mocouples eliminated silica based materials as gradient 
zone insulation candidates. 

The heat extraction plate material and thick- 
ness comprised the second option for maximizing the 
AADSF gradient capability for the USMP-2 experi- 
ment Several factors influenced the heat extraction 
plate design. Since the heat extraction plate conducts 
heat from the hot zone to the cooling system, the plate 
must withstand the thermal stresses induced by large 
radial temperature gradients. Radial slots were incorpo- 
rated into the plate design to reduce thermal stresses in- 
duced by radial temperature gradients. The need to 
accommodate axial thermal expansion of the cold zone 
heater core limited the allowable flexural stiffness of the 
heat extraction plate. This low stiffness requirement led 
to a multi-layered plate design, which provides the nec- 
essary flexibility without sacrificing the radial thermal 
conductance needed to produce high gradients. 

Both upper and lower bounds exist on the heat 
extraction plate conductance. Gradient performance re- 
quirements define the lower bound, and booster heater 
power limits define the upper thermal conductance 
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bound. Large plate conductances could degrade hot 
zone isothermality by loading the booster heater beyond 
its power limits. The relationship between the heat ex- 
traction plate and booster heater limits depends on the 
thermal conductance between the plate and booster 
heater, thus coupling the heat extraction plate and gra- 
dient zone insulation configurations to each other. 

Thermal Modeling and Analysis 

During development of the AADSF, several 
thermal models of the furnace were developed and used 
as design tools. The majority of the design supporting 
analysis was performed with an axi symmetric 
steady-state model of the furnace. This model was used 
to evaluate gradient zone configurations to ensure con- 
formance to design specifications and mission science 
goals. Also, the model provided estimated heater power 
requirements to support sizing of the heaters, and deter- 
mined component temperature distributions to support 
thermoelastic stress analysis. The model was also used 
to predict furnace performance and sample thermal pro- 
files for the USMP-2 crystal growth experiment. 

The AADSF thermal model is based on the fi- 
nite difference method, and was developed for solution 
with the Systems Improved Numeric Differencing Ana- 
lyzer (SINDA/G™) thermal solver. The S1NDA ther- 
mal resistance-capacitance network was developed with 
PATRAN™. Although PATRAN™ is primarily a fi- 
nite element mesh generator, the finite difference SIN- 
DA model was created by assigning lumped mass 
thermal nodes to each of the finite element centroids. 

Thermal modeling goals drove the choice of 
analysis methods and tools. The high degree of flexibil- 
ity that can be built into a SIND A model led to its selec- 
tion as the thermal solver. The need to perform 
parametric design studies necessitated a flexible ther- 
mal model permitting efficient specification of material 
properties (heat extraction plate and gradient zone in- 
sulation) and some component dimensions (heat extrac- 
tion plate thickness), without repetitive preprocessing. 
In addition to meeting these needs, the PC based ver- 
sion of SINDA provided an additional benefit. Once 
compiled, PC-SINDA models are portable and may be 
executed on any PC having a compatible processor. The 
resulting model can be made available to the Principal 
Investigator, and permits running simultaneous simu- 
lations on multiple PCs. 

A computer aided design (CAD) program was 
used to develop an idealized cross sectional diagram of 
the furnace. From preliminary thermal analysis, the 
gradient zone thickness for the USMP-2 experiment 
was selected to be 14 mm, with a 1.5 mm heat extrac- 
tion plate. This gradient zone configuration provided 


the geometric basis for the AADSF thermal model. 
Each furnace component was represented, along with 
the cartridge, ampoule, and sample material. The car- 
tridge was positioned in the AADSF with the sample 
approximately centered in the gradient zone. The ther- 
mal model geometry corresponds to this assumed car- 
tridge position. 

Using the CAD program, each region repre- 
senting a component or gas filled cavity was subdivided 
to better control meshing operations. When resistance- 
capacitance networks are defined by computing the 
thermal resistances between finite element centroids, lo- 
cal numerical errors result in the neighborhood of un- 
structured or nonrectangular mesh elements. Additional 
CAD based premeshing operations minimized the num- 
ber of nonrectangular elements, and ensured that un- 
avoidable nonrectangular elements occurred only in 
noncritical areas of the model. The thermal model ge- 
ometry and computational mesh are shown in Figure 4. 



Geometric Idealization of Axtsym metric Two Dimensional 

the AADSF Computational Mesh 

Figure 4: AADSF Thermal Model Geometry and Mesh 

The CAD description of model geometry was 
imported into PATRAN™, which was then used to gen- 
erate the detailed computational mesh. This mesh con- 
sists principally of rectangular elements, but includes 
some nonrectangular quadrilateral, triangular, and line 
elements. PATRAN™ was also used to generate 
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boundary conditions and material property assignments. 
The resulting finite element mesh was translated to a 
resistance-capacitance network in SIND A format. The 
SINDA input file resulting from this operation con- 
tained a complete description of the thermal conduction 
portion of the problem, but lacked material property 
values and radiation heat transfer modeling for cavities 
internal to the furnace. 

Radiation conductance networks were com- 
puted for all gas filled voids inside the furnace. Such 
cavities include the annular regions separating the heat- 
er cores and primary insulation, the cartridge and fur- 
nace bore, and the sample/ampoule and cartridge tube. 
During mesh development, PATRAN™ was used to 
generate line elements on the surfaces of these regions. 
These elements provided a convenient geometric basis 
from which to construct thermal radiation models. En- 
gineering judgment and a literature review [3, 4, 5) led 
to best estimates for the total hemispherical emissivity 
of the various radiating surfaces at furnace operational 
temperatures. These emissivity data completed the mod- 
el description for computation of the thermal radiation 
conductance networks. 

The Net Energy Verification and Determina- 
tion Analyzer (NEVADA™) program was used to com- 
pute diffuse gray body solutions for the radiation 
conductance networks. View factors were statistically 
integrated using the Monte Carlo method. A total of 
50,000 rays were randomly generated from each sur- 
face. The view factors were numerically integrated by 
counting the intersections of these rays with other sur- 
faces in the model. With 50,000 rays, view factors as 
small as 0.01 are accurate to ±10% with 95% confi- 
dence [6], The error bands decrease for larger view fac- 
tors. Once the view factors were computed, gray body 
solutions were obtained with the Gebhart solver [7] in- 
cluded with NEVADA™ and translated to SINDA 
format. 

Thermal conductivity data remained to be add- 
ed to the AADSF thermal model. A literature review [3, 
4, 5, 8, 9] resulted in arrays of temperature dependent 
thermal conductivity data. Addition of the radiation 
solutions and material data completed the AADSF ther- 
mal model description. The resulting model contained 
2,150 nodes and 12,166 conductances interconnecting 
these nodes. Radiation modeling accounted for 68% of 
the conductors. 

Although the AADSF thermal model descrip- 
tion was essentially complete at this point, it was not 
ready for use as a parametric analysis tool. A propor- 
tional integrating control algorithm was added to the 
model to automatically adjust power to achieve the de- 
sired heater set point temperatures during iterative 
steady state solutions. On each iteration, the controller 


computes power adjustments from the set point error 
and the change in error since the previous iteration. 
Code was also added to permit specification of set point 
temperatures, material properties, and heat extraction 
plate thickness without requiring repetitive preprocess- 
ing. With these enhancements, the AADSF thermal 
model became a useful tool for parametric analysis. 

Two major limitations apply to the AADSF 
thermal model. First, the model simulates only steady 
state conditions. Transient analysis was not needed to 
support the design effort. Second, since the AADSF 
thermal model is based on the differential equation of 
heat conduction, only results derivable from tempera- 
tures and heat fluxes can be computed. Convection and 
diffusion effects within the liquid portion of the sample 
cannot be simulated with this model. 

Thermal Analysis and Application to Design 

The AADSF thermal model served two main 
purposes. First it served as a parametric design analysis 
tool. The model was used to characterize AADSF gradi- 
ent performance as a function of operating temperature 
and gradient zone design variables, including gradient 
zone insulation material and heat extraction plate con- 
figuration. The calculated peak axial temperature gradi- 
ent induced in a quartz rod provided an analytical 
measure of gradient performance. The furnace model 
was also used to predict AADSF performance for the 
USMP-2 experiment. 

Over 270 simulations were performed for the 
parametric analysis. For all cases, the thermal conduc- 
tivity of quartz was used to represent the sample. In ef- 
fect, this modeled a uniform quartz rod within the 
metallic cartridge, and provided a means of comparing 
thermal model predictions with performance goals and 
test data. Five parameters were varied in the study: (1) 
hot zone operating temperature - from 400°C to 
1200°C, (2) cold zone control mode - isothermal or 
floating, (3) heat extraction plate material - Inconel 718 
or pure platinum, (4) heat extraction plate thickness - 
ranging from zero to 3 mm, and (5) gradient zone in- 
sulation material - dense alumina, AL-30™, or FBD™. 

When considering high thermal gradient op- 
erations, two bounding cold zone operational conditions 
are of interest These two conditions are: (1) isothermal 
cold zone - the cold main and cold guard heaters are 
commanded to the same temperature, and (2) floating 
cold zone - the cold zone idles at minimum power with 
the cold guard heater temperature below that of the cold 
main isothermal cold zone operation requires signifi- 
cant power to the cold guard heater. Floating cold zone 
operation requires more heat from the booster heater for 
any given hot zone temperature, but typically results in 
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slightly less overall furnace power consumption. For 
true floating cold zone operation, the cold main and 
cold guard heaters would consume zero power. Howev- 
er a minimal power level, equivalent to one ampere to 
each heater, was applied for this analysis. This provided 
better estimates of minimum controllable cold zone 
temperatures for non-isothermal operation, and had a 
conservative effect on gradient performance predictions. 
Analysis clearly showed that, at the expense of cold 
zone isothermality, floating cold zone operation signifi- 
cantly increases the gradient capability of the AADSF. 

For each steady state simulation, a converged 
solution was achieved using a successive over- 
relaxation method [10] to solve for all nodal tempera- 
tures. The solutions were not considered converged un- 
less two criteria were met: (1) all nodal temperatures 
vary by less than 0.00 1°C between successive iterations, 
and (2) the net system energy imbalance was less than 
1%. For each of 270 cases, the nodal temperature pro- 
file was computed. Results of specific interest, such as 
heater power requirements, cold zone temperatures, and 
peak axial quartz rod temperature gradients were 
computed. 

The AADSF thermal model was run for five 
different hot zone temperatures ranging from 400°C to 
1200°C in 200°C increments. An isothermal hot zone 
was simulated by commanding the hot guard, hot main, 
and booster heaters to the same set point temperature, 
and the peak axial temperature gradient induced in a 
quartz rod was computed. In general, larger gradients 
were produced at higher hot zone temperatures. The 
greater temperature difference between the hot zone and 
the cooling system allows larger temperature drops 
across the gradient zone to be developed, and thus larg- 
er gradients are produced. Gradient performance is fur- 
ther augmented at higher temperatures because of 
enhanced radiation heat transfer. Higher temperatures 
produce greater thermal radiation conductance, which 
results in sample temperature variations that correspond 
more directly with the steep thermal gradients of the 
furnace bore. 

Both isothermal and floating cold zone opera- 
tional modes were considered in this parametric analy- 
sis. Figure 5 plots the predicted axial gradient induced 
in a quartz rod versus hot zone temperature for a 1.5 
mm Inconel 718 heat extraction plate and AL-30™ gra- 
dient zone insulation. The thermal model indicates that 
floating cold zone operation can more than double the 
AADSF gradient capability. For hot zone operation at 
1100°C, a gradient in excess of 125°C/cm is predicted. 
Isothermal operation produces only 61°C/cm, short of 
the 100°C/cm goal. These results indicate that the first 
flight gradient objectives require operating the AADSF 



Figure 5: AADSF Gradient Performance versus Coid Zone 
Control Mode 


with the cold guard heater commanded to a lower 
temperature than the cold main heater. 

The parametric analysis considered two candi- 
date high temperature materials for the heat extraction 
plate: platinum and Inconel 718. Figure 6 plots the pre- 
dicted axial gradient induced in a quartz rod versus the 
hot zone temperature for a 1.5 mm heat extraction plate 
and AL-30™ gradient zone insulation, for both con- 
trolled and floating cold zone operation. For isothermal 
cold zone operation, a platinum plate significantly 



Figure 6 : AADSF Gradient Performance versus Heat Ex- 
traction Plate Material and Cold Zone Control Mode 
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improves gradient performance, but still fails to meet 
the 100°C/cm required gradient for a hot zone at 
1100°C. For floating cold zone operation, the platinum 
heat extraction plate results in somewhat higher gradi- 
ents than an Inconel 718 plate, although both materials 
produce gradients in excess of the 100°C/cm design 
goal. As the hot zone temperature increases, the gradi- 
ent performance gains that can be achieved with a plat- 
inum plate become more pronounced for both cold zone 
control modes. Since the USMP-2 experiment will be 
operated with the cold guard heater temperature below 
that of the cold main. Inconel 718 was selected for the 
heat extraction plate, because it met the gradient goal at 
a lower cost. Future AADSF experiments, if intended 
for high gradient and high temperature operation 
(above 1150°C), or if cold zone isothermality is more 
critical, may employ platinum heat extraction plates. 

A range of candidate thicknesses were consid- 
ered for the heat extraction plate in the design analysis: 
from zero to 3 mm in 0.75 mm increments. Figure 7 
plots the predicted axial gradient induced in a quartz 



Figure 7; AADSF Gradient Performance versus Heat 
Extraction Plate Thickness for Isothermal Cold Zone 
Operational Mode 


rod versus the hot zone temperature for isothermal cold 
zone operation, for the furnace equipped with an In- 
conel 718 heat extraction plate and AL-30™ gradient 
zone insulation. Although the 100°C/cm goal is not 
met, these results show that the heat extraction plate 
thickness can have a profound influence on gradient 
performance. Figure 8 plots the same results for floating 
cold zone operation. These results indicate that the 
100°C/cm gradient goal is met for a hot zone operating 
at 1 100°C without a heat extraction plate. The thermal 



Hot Zone Temperature (*C] 

Figure 8: AADSF Gradient Performance versus Heat 
Extraction Plate Thickness for Floating Cold Zone 
Operational Mode 

model shows that for floating cold zone operation, gra- 
dient performance is somewhat less sensitive to the heat 
extraction plate thicknesses than for isothermal cold 
zone operation. During floating cold zone operation, the 
cold guard heater operates at a lower temperature than 
the cold main heater. The temperature drop between 
these heaters provides an effective means of conducting 
hot zone heat transferred across the gradient zone to the 
cold zone end support structure. In the presence of this 
thermal path, the heat extraction plate thickness be- 
comes less important with respect to gradient perform- 
ance. However, at higher operating temperatures, the 
thermal conductivity simultaneously increases for the 
metallic heat extraction plate and decreases for the ber- 
yllia heater core materials, so an increase in heat ex- 
traction plate thickness produces more pronounced 
gradient performance gains. 

Three candidate materials were considered for 
the gradient zone ins ulation: AL-30™, FBD™, and 
dense alumina. AL-30™ exhibits the best insulating 
properties, and dense alumina has the greatest mechani- 
cal strength. Figure 9 plots the predicted axial gradient 
induced in a quartz rod versus hot zone temperature for 
a 1.5 mm Inconel 718 heat extraction plate and floating 
cold zone operation. For a hot zone at 1100°C, 
AL-30™ and FBD™ both meet the 100°C/cm gradient 
goal, while dense alumina produces a gradient of only 
52°C/cm. AL-30™ gradient zone insulation produces 
only slightly larger gradients than FBD™, despite its 
significantly lower thermal conductivity. This is due to 
heat conduction through the metallic insert that me- 
chanically connects the gradient zone and heater cores 
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Hot Zone Temperature [*CJ 

Figure 9: AADSF Gradient Performance versus Gradient 
Zone Insulation Material 


to the furnace structure through the heat extraction 
plate. This insert provides a conduction path parallel to 
the gradient zone insulation, and renders the AADSF 
gradient performance relatively insensitive to the gradi- 
ent zone insulation thermal conductivity values compa- 
rable to AL-30™ and FBD™. AL-30™ was chosen 
over the FBD™ material primarily because FBD™ 
principally contains zirconia, which becomes electrical- 
ly conductive at high temperatures. This electrical con- 
ductivity could produce unwanted coupling between the 
resistive heater elements and their control 
thermocouples. 

Performance Predictions for First Flight 
Confieuration 

Once the gradient zone was defined with 
AL-30™ insulation and a 1.5 mm Inconel 718 heat ex- 
traction plate, the AADSF thermal model was used to 
predict furnace performance for the USMP-2 experi- 
ment Three cases were considered: (1) maximum hot 
zone and cold zone temperature, (2) maximum hot zone 
temperature and maximum gradient, and (3) USMP-2 
experiment conditions. 

The first case simulated maximum temperature 
operation in all heated zones with the hot zone at 
1150°C and the cold zone at 850°C. The second case 
simultaneously simulated 1150°C maximum hot zone 
temperature and maximum gradient conditions. These 
two cases were analyzed to check heater power require- 
ments against design heat dissipation limits. For both 


cases, the cartridge was assumed to contain a quartz rod 
rather than a sample/ampoule assembly. The cold main 
heater was predicted to float over its set point tempera- 
ture for both cases. Set point errors of 1.6°C and 41.0°C 
were predicted for 850°C and 600°C cold zones, respec- 
tively. For the second case, a peak axial temperature 
gradient of 105°C/cm was predicted for the quartz rod. 
Table 2 lists the predicted heater power requirements. 
No power exceedances resulted from comparing the 
heat dissipated within the cores to the design limits. 
The results also indicate hot guard and hot main power 
requirements to be largely independent of cold zone 
temperature. However, reducing the cold zone set point 
from 850°C to 600°C significantly increased booster 
heater power demand. 


! Design 

Thermal Analysis Results 

Heater 

Desig- 

nation 

Core 

Power 

Limit 

(Watts) 

115Q/850°CHot/CoW 

Zone 

1150/600°C Hot/Coki 
Zone 

Core 

Power 

(Watts) 

Total 

Power 

(Watts) 

Core 

Power 

(Watts) 

Total 

Power 

(Watts) 

Hot Guard 

320 

228.5 

249.8 

227.4 

248.6 

Hot Main 

930 

159.9 

173 

164.9 

178.4 

Booster 

275 

117.5 

164.2 

149.1 

208.4 

Cold Main 

130 

0 

0 

0 

0 

Cold 

Guard 

195 

1653 

188.9 

64.1 

74 

Total 

671.2 

775.9 

605.5 

709.4 


Table 2: Predicted AADSF Power Performance at Design 
Temperature Limits 


The third analysis case simulated AADSF per- 
formance for the USMP-2 experiment. For this simu- 
lation, all three hot zone heaters were commanded to 
850°C, the cold main heater to 350°C, and the cold 
guard heater to 325°C. Table 3 lists the predicted heater 
power requirements for this operational condition. The 
results indicate the furnace should operate within its de- 
sign power limits for the USMP-2 experiment. All five 
heaters were predicted to be powered and to operate 
without set point errors, indicating that active thermal 
control to these temperatures is possible for all five 
heated zones. Figure 10 illustrates the predicted axial 
temperature and temperature gradient profile through 
the sample. For this analysis, a representative melting 
temperature of 700°C was assumed. A sudden sample 
thermal conductivity change at the solidification tem- 
perature provided a simple means of representing the 
phase change at the interface. Actually, the thermo- 
physical behavior of the USMP-2 sample material is 
considerably more complex. During processing, the 
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Design 

Thermal Analysis Results 

Heater 

Designation 

Core Power 
Limit (Watts) 

850°C Hot Zone and 35G/325°C 
Cold Main/Guard 

Core Power 
(Watts) 

Total Power 
(Watts) 

Hot Guard 

320 

179.4 

197 

Hot Main 

930 

98.9 

107.3 

Booster 

275 

104.6 

146.8 

Cold Main 

130 

7.2 

8.5 

Cold Guard 

195 

8.7 

10.2 

Total 

398.8 

469.8 


Table 3: Predicted USMP-2 AADSF Steady Stale Power 
Performance 



representation of the thermal paths throughout the fur- 
nace would be indicated by accurate total and individual 
heater power predictions. With reliable power predic- 
tions, the model can be used to evaluate the feasibility 
of various experimental scenarios. Heater power can be 
evaluated against design power limits, and uncontrolla- 
ble set point temperature combinations can be avoided. 

Figures 11 and 12 compare predicted and mea- 
sured steady state heater power for maximum power 
and maximum gradient operation, respectively. The re- 
sults in these figures do not quite match the data in Ta- 
ble 2, because a customized empty-bore version of the 
model was used to better match the test conditions. 
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Figure 11: Predicted versus Measured Power for Maximum 
Hot Zone and Cold Zone Temperature 


Figure 10: Predicted Sample Axial Temperature and Ther- 
mal Gradient Profile for the USMP-2 Experiment 

sample melt becomes enriched with solute near the in- 
terface, resulting in compositional variations that cause 
the solidification temperature to change. The results in 
Figure 10 show a peak axial gradient in the sample of 
106°C/cm. The location of the peak gradient corre- 
sponds closely to the 700°C melting point in the 
sample. 

Comparison to Test 

Steady state heater power predictions obtained 
from the AADSF thermal model have been compared to 
data obtained during ground testing. Power predictions 
provide the most effective means of assessing the pre- 
dictive capabilities of this model, for two reasons. First, 
an understanding of power requirements is fundamental 
to the furnace design process. Second, a reliable 



11 SOX 11 SOX 1150X 600 X 600 X 


Figure 12: Predicted versus Measured Power for Maximum 
Hot Zone Temperature and Maximum Gradient 
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Predicted and measured distributions are similar, with 
predicted total power less than 5% higher than mea- 
sured values. For the high gradient condition, the model 
predicted the cold main heater to operate 49°C above its 
600°C set point. In the test, the steady state cold main 
temperature was 640°C. 

Figure 13 compares predicted and measured 
power data for conditions representative of the USMP-2 
experiment. The predicted and measured power distri- 
bution are quite similar, with predicted total power ex- 
ceeding measured power by 8%. In Figure 13, the 
model more accurately predicts the power distribution 
than the empty-bore version shown in Figures 1 1 and 
12. This is attributed to natural convection in the fur- 
nace bore during ground testing. The presence of the 
cartridge suppresses this convection during testing. 
Since the model includes no bore convection effects, 
predictions should correlate more closely to measure- 
ments for cartridge- in-bo re conditions than for empty- 
bore conditions. 



Figure 13: Predicted versus Measured Power for 
First Flight Conditions 

Figure 14 compares measured and predicted 
temperature profiles within the cartridge. The test was 
performed with an alumina rod instrumented with four 
thermocouples, mounted inside a cartridge. Alumina 
thermal conductivity values were assigned to the region 
representing the sample in the model to better simulate 
this test condition. Figure 14 shows two presentations of 
the test data: (1) the thermocouple temperature closest 
to the cold end of the alumina rod recorded throughout 
cartridge translation, and (2) four individual tempera- 
ture measurements taken at the end of translation. The 
temperature distribution predicted with the thermal 
model agrees well with measurements, although the 


model conservatively under-predicts the temperature 
gradient. 



Figure 14: Predicted versus Measured Temperature for 
Calibration Sample and USMP-2 Temperatures 

Conclusions 

Thermal modeling proved to be an effective 
design tool during development of the AADSF. A ther- 
mal model was specifically developed to perform para- 
metric analysis and used to characterize how the 
gradient zone design affects the thermal performance of 
the furnace. The results of the parametric analysis were 
used to select a gradient zone configuration to fulfill the 
requirements of the USMP-2 experiment. Furthermore, 
furnace characterization was achieved in far less time 
and expense with the thermal model than would be re- 
quired to build a prototype and run an equivalent series 
of characterization tests. 

Design analysis of the AADSF system predict- 
ing successful USMP-2 performance was later verified 
during tests conducted at Marshall Space Flight Center. 
Overall, the results obtained from the AADSF model 
compared well to measurements taken during tests. 
Steady state heater power requirements and temperature 
measurements obtained with a quartz rod instrumented 
with thermocouples were often within a few percent of 
predictions obtained with the thermal model. The ap- 
parent accuracy of the model is attributed largely to 
thorough handling of thermal radiation and temperature 
dependent material properties. 

The high degree of flexibility built into the 
AADSF thermal model resulted in its success as a 
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parametric analysis and design tool. Code was added to 
permit user specification of design parameters such as 
the temperature dependent component thermal conduc- 
tivity data or the heat extraction plate thickness. More 
importantly, code was added to automatically adjust 
heater powers to achieve the specified set point tem- 
peratures. Without this numerical heater control, heater 
powers would have to be specified and adjusted based 
on resulting temperatures, a process that would have 
rendered parametric design analysis impractical. 

Several future expansions and uses for the 
AADSF thermal model are suggested by this work. 
With a thorough set of furnace characterization test 
data, the thermal model could be fine tuned to enhance 
its predictive capabilities. The model could be upgraded 
to transient capability and used to develop experiment 
timelines and fine tune heater control algorithms. The 
model could also aid in detailed analysis of crystal 
growth within the ampoule by providing realistic ther- 
mal boundary conditions at the furnace bore or car- 
tridge wall. Because the AADSF thermal model is 
portable, it is available to the investigator to help devel- 
op effective crystal growth experiments. The model can 
be used to assess a large number of processing profiles 
before performing actual experiments, thereby saving 
valuable laboratory resources and experiment develop- 
ment time. 
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